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Abstract 

The existence of anticipatory, complete and lag synchronization in a single system having two 
different time-delays, that is feedback delay t\ and coupling delay T2, is identified. The transition 
from anticipatory to complete synchronization and from complete to lag synchronization as a 
function of coupling delay t<i with suitable stability condition is discussed. In particular, it is 
shown that the stability condition is independent of the delay times t\ and t<i- Consequently for 
a fixed set of parameters, all the three types of synchronizations can be realized. Further the 
emergence of exact anticipatory/complete/lag synchronization from the desynchronized state via 
approximate synchronization, when one of the system parameters (62) is varied, is characterized 
by the minimum of similarity function and the transition from on-off intermittency via periodic 
structure in laminar phase distribution. 
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I. INTRODUCTION 



Synchronization phenomena dates back to the period of Huygens in 1665, when he found 
that two very weakly coupled pendulum clocks, hanging from the same beam, become 
phase synchronized l|. Since the early identification of synchronization in chaotic oscil- 

QDD 



lators 
areas of science 



;he phenomenon has attracted considerable research activity in different 
irfa and Severn, generations and testing appneations have been 
developed. Chaos synchronization phenomenon is of interest not only from a theoretical 
point of view but also has potential applications in diverse subjects such as as biological, 
neurological, laser, chemical, electrical and fluid mechanical systems as well as in secure 

communication, cryptography and so on. DDDDDODDQQ 



A recent review 

on the phenomenon of chaos synchronization can be found in the reference jllj]. In recent 
years, different kinds of synchronization have been identified: Complete (or identical) syn- 
chronization gene ralized synchronization 

lag synchronization llg, [19J and anticipatory synchronization [2(, 



13l Il4j , phase synchronization 
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For a critical 

discussion on the interrelationship between various kinds of synchronization, we may refer 
to refs. Transition from one kind of synchronization to the other, coexistence of 

different kinds of synchronization in time series and also the nature of transition have also 



been studied extensively 







18, 



19 



25 



'), 26] in coupled chaotic systems. 
One of the most important applications of chaos synchronization is secure communica- 
tion. It is now an accepted fact that secure communication based on simple low dimensional 
chaotic systems does not ensure sufficient level of security, as the associated chaotic attrac- 
tors can be reconstructed with some effort and the hidden message can be retrieved by an 
eavesdropper 2jj. One way to overcome this problem is to consider chaos synchronization 
in high dimensional systems having multiple positive Lyapunov exponents. This increases 
security by giving rise to much more complex time series, which are apparently not vul- 
nerable to the unmasking procedures generally. Recently chaotic time-delay systems have 
been suggested as good candidates for secure communication [23, Q, as the time-delay 
systems are essentially infinite dimensional in nature and are described by delay differential 
equations, and that they can admit hyperchaotic attractors with large number of positive 
Lyapunov exponents for suitable nonlinearity. Therefore the study of chaos synchronization 
in the time-delay systems is of considerable practical significance. However, it should be 
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noted that one has to be cautious due to the fact that even in time-delay systems with 
multiple positive Lyapunov exponents unmasking may be possible. Particularly, this is so 
if any reconstruction of the dynamics of the system is achieved in some ap prop riate space 
even for very high dimensional dynamics as demonstrated by Zhou and Lai [30[ in the case 
of Mackey- Glass equation. 

Time-delay is ubiquitous in many physical systems due to finite switching speed of am- 
plifiers, finite signal propagation time in biological networks, finite chemical reaction times, 



28, 



29, 31, 32]. In recent times, considerable work has been 



memory effects and so on 

carried out on the effect of time-delay in limit cycle oscillators [3, time-delay feed- 
back j2, 0], networks with time-delay coupling j^, etc. Recently, we have shown that 
even a single scalar delay equation with piecewise linear function can exhibit hyperchaotic 



behavior even for small values of time-delay 



It is therefore of importance to consider 



the synchronization of chaos in such scalar piecewise linear delay differential systems with 
appropriate delay coupling. Interestingly, in the present work we find that in such a coupled 
system, one can identify anticipatory, complete and lag synchronizations by simply tuning 
the second time-delay parameter in the coupling, for a fixed set of system parameters sat- 
isfying appropriate stability condition. The results have been corroborated by the nature 
of similarity functions, and transition behavior characterized by probability distribution of 
laminar phase during approximate synchronization which precedes the exact synchronization 
when a system parameter is varied. We also wish to point out that to our knowledge such 
transitions between all the three types of above synchronizations have not been reported in 
nonhyperchaotic systems and it appears that the present type of hyperchaotic systems are 
convenient tools to realize such transitions by tuning the delay parameters suitably. 

Specifically, in this paper we will consider chaos synchronization of two single scalar 
piecewise-linear time-delay systems studied in the references j^l) 38, 3^| with unidirectional 
coupling between them and having two different time-delays: one in the coupling term and 
the other in the individual systems, namely, feedback delay. We have identified the stability 
condition for synchronization following Krasovskii-Lyapunov theory and demonstrate that 
there exists transition between three different kinds of synchronization, namely anticipatory, 
complete and lag synchronizations, as a function of time-delay in the coupling. To charac- 
terize the existence of anticipatory and lag synchronizations, we have plotted the similarity 
function S(t). We have also demonstrated that when the system parameter 62 is varied, 
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the onset of exact anticipatory/complete/lag synchronization from the desynchronized state 
is preceded by a region of approximate synchronized state. We also show that the later is 
characterized by a transition from on-off intermittency to aperiodic structure in the laminar 
phase distribution, as suggested in the work of Zhan et al. |l9( for the case of lag synchroniza- 
tion. The plan of the paper is as follows. In sec. II, we introduce the unidirectionally coupled 
scalar time-delay system and identify the condition for stability of synchronized states. In 
sec. Ill, we point out the existence of anticipatory synchronization when the strength of 
the coupling delay is less than feedback delay, while in sec. IV, complete synchronization is 
realized when the two delays are equal. Lag synchronization is shown to set in when the 
coupling delay exceeds the feedback delay in sec. V. Finally in sec. VI, we summarize our 
results. 



II. PIECEWISE LINEAR TIME-DELAY SYSTEM AND STABILITY CONDI- 
TION FOR CHAOS SYNCHRONIZATION 

At first, we will introduce the single scalar time-delay system with piecewise linearity 
and bring out its hyperchaotic nature for suitable values of the system parameters. Then 
the unidirectional delay coupling is introduced between two scalar systems and the stability 
condition for chaos synchronization is derived. 



A. The scalar delay system 



We consider the following first order delay differential equation introduced by Lu and 



He 



31] 



and discussed in detail by Thangavel et al. 

x{t) = -ax(t) +bf(x(t-r)), (1) 

where a and b are parameters, r is the time-delay and / is an odd piecewise linear function 
defined as 

0, x < -4/3 

-1.5a; -2, -4/3 < x < -0.8 

x, -0.8 < x < 0.8 (2) 

-1.5x + 2, 0.8<x<4/3 

0, x > 4/3 



/(*) 



It is also of interest to consider additional forcing on the right hand side of Eq. (JI}; however, 
this is not considered here. The schematic form of © is shown in Fig. 1. Recently, we 
that systems of the form exhibit hyperchaotic behavior for suitable 
parametric values. For our present study, we find that for the choice of the parameters a = 
1.0,6 = 1.2 and r = 25.0 with the initial condition x(t) = 0.9, t e (-5,0), Eq. (JTJ) exhibits 
hyperchaos. The corresponding pseudoattractor is shown in the Fig. 2. The hyperchaotic 
nature of Eq. is confirmed by the existence of multiple positive Lyapunov exponents. 
The first ten maximal Lyapunov exponents for the parameters a = 1.0,6 = 1.2, x(t) = 
0.9, t G (—5,0), as a function of time-delay r is shown in Fig. 3, which are evaluated using 
the procedure suggested by J. D. Farmer [29]. 



B. Coupled system and the stability condition 

Now let us consider the following unidirectionally coupled drive x\(t) and response X2 it) 
systems with two different time-delays Tl and r 2 as feedback and coupling time-delays, 
respectively, 

x\(t) = -axi(t) + b 1 f{x 1 {t - n)), (3a) 
x 2 (t) = -ax 2 (t) + b 2 f(x 2 (t - n)) + fea/fai (< - r 2 )), (3b) 

where 61, 6 2 and 63 are constants, a > 0, and f(x) is of the same form as in Eq. (J2J). 

Now we can deduce the stability condition for synchronization of the two time-delay 
systems Eqs. (JH^) and (JUb) in the presence of the delay coupling 6 3 /(xi(t — r 2 )). The time 
evolution of the difference system with the state variable A = x\ T2 - Tl — x 2 (which corresponds 
to anticipatory synchronization when r 2 < Ti, identical synchronization for t 2 = T\ and lag 
synchronization when r 2 > Ti), where Xi T2 „ T1 = X\(t — (r 2 — Ti)), can be written for small 
values of A by using the evolution Eqs. (3) as 

A = -aA + (6 2 + 63 - b l )f(x 1 (t - r 2 )) + b 2 f{x l {t - r 2 ))A Tl , A T = A(t - r) (4) 

In order to study the stability of the synchronization manifold, we choose the parametric 
condition, 

61 = 6 2 + 63, (5) 
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so that the evolution equation for the difference system A becomes 

A = -aA + b 2 f'(x 1 (t-r 2 ))A Tl . (6) 

The synchronization manifold is locally attracting if the origin of this equation is stable. 
Following Krasovskii-Lyapunov functional approach [39,40], we define a positive definite 
Lyapunov functional of the form 

l r° 

V(t) = -A 2 + /2 / A 2 (t + 6)d6, (7) 

* J-Tl 

where \x is an arbitrary positive parameter, \i > 0. Note that V(t) approaches zero as A — > 0. 

To estimate a sufficient condition for the stability of the solution A = 0, we require the 
derivative of the functional V(t) along the trajectory of Eq. (JHJ), 

^ = -aA 2 + b 2 f( Xl (t - r 2 ))AA Tl + M A 2 - fiA 2 n , (8) 

to be negative. The above equation can be rewritten as 

^ = -/iA 2 r(X, yU ), (9) 



where X = A n /A, T =[((a - fi)/fi) -{b 2 f{xi{t - t 2 ))/ll)X + X 2 ]. In order to show that 
^ < for all A and A r and so for all X, it is sufficient to show that Y min > 0. One 
can easily check that the absolute minimum of T occurs at X = ^hf'ix^t - r 2 )) with 
Tmin =[4/i(a — //) — b\f\x\{t — t 2 )) 2 ] /4/x 2 . Consequently, we have the condition for stability 

as 

a>^f'(x 1 (t-T 2 )) 2 + Li = ^). (10) 

Again as a function of /x for a given /'(x) has an absolute minimum at fi — (|& 2 /'(a;i(t — 
r 2 ))|)/2 with $ min = \b 2 f'(xx(t - r 2 ))|. Since $ > $ min = |6 2 /' ( Xl (t - r 2 ))|, from the 
inequality (fTHjl . it turns out that the sufficient condition for asymptotic stability is 

a>|6 2 /'(a; 1 (t-r 2 ))| (11) 

along with the condition (jSJ) on the parameters b\, b 2 and 63. 

Now from the form of the piecewise linear function f(x) given by Eq. (2), we have, 

. 1.5, 0.8 < < I 
\f\x 1 {t-r 2 ))\ = { ~ 11 " 3 (12) 

1.0, xi < 0.8 
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Note that the region \x±\ > 4/3 is outside the dynamics of the present system (see Eq. ©). 
Consequently the stability condition (jllj) becomes a > 1 - 5 1 £»2 1 > I&2I along with the para- 
metric restriction bi = b 2 + 63. 

Thus one can take a > \b 2 \ as a less stringent condition for to be valid, while 



as the most general condition specified by (jllj) for asymptotic stability of the synchronized 
state A = 0. The condition (fT3*J) indeed corresponds to the stability condition for exact 
anticipatory, identical as well as lag synchronizations for suitable values of the coupling 
delay r 2 . It may also be noted that the stability condition (|T3~j) is independent of the both 
the delay parameters t% and t 2 . In the following, we will demonstrate the transition from 
anticipatory to lag synchronization via complete synchronization as the coupling delay r 2 is 
varied from t 2 < t\ to r% > tl, subject to the stability condition (fT3*J) with the parametric 
restriction b\ = b 2 + 63. However, we also point out from detailed numerical analysis that 
when the less general condition X .5 1 &2 1 > a > I b 2 | is satisfied, approximate synchronization 
(anticipatory/complete/lag) occurs. 

III. ANTICIPATORY SYNCHRONIZATION FOR t 2 < n 

To start with, we first consider the transition to anticipatory synchronization in the 
coupled system Q. We have fixed the value of the feedback time-delay T\ at T\ = 25.0 while 
the other parameters are fixed as a = 0.16, 61 = 0.2, b 2 = 0.1,63 = 0.1 and the time-delay 
in the coupling r 2 is treated as the control parameter. With the above mentioned stability 
condition ()13|) and with the coupling delay t 2 being less than the feedback delay T\, one can 
observe the transition to anticipatory synchronization. The time trajectory plot is shown 
in Fig. 4a depicting anticipatory synchronization, for the specific value of r 2 = 20.0 with 
the anticipating time equal to that of difference between feedback and the coupling delays, 
that is, r = r 2 — t\. The time-shifted plot Fig. 4b, x 2 (t — r) Vs Xi(t), shows a concentrated 
diagonal line confirming the existence of anticipatory synchronization (We may note here 
that in all our numerical studies in this paper we leave out sufficiently large number of 
transients, before presenting our figures). 



Sometime ago, Rosenbulm et al. |15( have introduced the notion of similarity function 



a > 1.516. 



'2|, 



(13) 
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Si(r) for characterizing the lag synchronization as a time averaged difference between the 
variables x\ and x 2 (with mean values being subtracted) taken with the time shift r, 



where, (x) means time average over the variable x. If the signals x\(t) and X2(t) are indepen- 
dent, the difference between them is of the same order as the signals themselves. If X\(t) = 
X2(t), as in the case of complete synchronization, the similarity function reaches a minimum 
S(t) = for r = 0. But for the case of nonzero value of time shift r, if Si(t) = 0, then there 
exists a time shift r between the two signals xi(t) and x^{t) such that x%{t + r) = xi(t), 
demonstrating lag synchronization. 

In the present study, we have used the same similarity function S/(t) to characterize 
anticipatory synchronization with negative time shift — r instead of the positive time shift 
r in Eq. (|14j) . In other words, one may define the similarity function for anticipatory 
synchronization as 



Then the minimum of S' a (r),that is S a (r) = 0, indicates that there exists a time shift 
—t between the two signals x\(t) and x 2 (t) such that x 2 (t ~ T ) = x i(t), demonstrating 
anticipatory synchronization. Fig. 5 shows the similarity function S a (r) as a function of 
the coupling delay T2 for four different values of 62, the parameter whose value determines 
the stability condition given by Eq. (|13p. while satisfying the parametric condition b\ = 
b 2 + &3- Curves 1 and 2 are plotted for the values of 62 = 0.18(> a = 0.16 > a/1.5) and 
62 = 0.16(= a > a/1.5), respectively, where the minimum values of S a (r) is found to be 
greater than zero, indicating that there is no exact time shift between the two signals x\(t) 
and x 2 (t). Note that in the both cases the stringent stability condition (fTHj) and the less 
stringent condition a > \b 2 \ are violated. Curve 3 corresponds to the value of b 2 = 0.15 
(which is less than a but greater than a/1.5), where the minimum value of S a (r) is almost 
zero, but not exactly zero (as may be seen in the inset of Fig. 5), indicating an approximate 
anticipatory synchronization X\(t) pa Xz(t — r). On the other hand the curve 4 is plotted 
for the value of 62 = 0.1(< a/1.5), satisfying the general stability criterion, Eq. (JT3J). It 
shows that the minimum of S a {r) = 0, thereby indicating that there exists an exact time 
shift between the two signals demonstrating anticipatory synchronization. The anticipating 




([x 2 (t + r) - Xl (t)} 2 } 

i(xi(t))(xm)] i/2 



(14) 




([x 2 (t - r) - x.m 

[(x\{t))(xl{t))\^ 



(15) 
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time is found to be equal to the difference between the coupling and feedback delay times, 
that is, r = r 2 — T\. Note that S a (r) = for all values of r 2 < n, indicating anticipatory- 
synchronization for a range of delay coupling. A further significance is that the anticipating 
time r = |r 2 — T\\ is an adjustable quantity as long as r 2 < Ti, which can be tuned suitably 
to satisfy experimental situations. 

Next, we show that the emergence of exact anticipatory synchronization is preceded by a 
region of approximate anticipatory synchronization, which is associated with the transition 
from on-off intermittency to a periodic structure in the laminar phase distribution as a 
function of the parameter 6 2 . First we choose the value of 6 2 as 6 2 = 0.17 (with b\ = 0.2 and 
63 = 0.03), above the value of a = 0.16, such that the general stability criterion, Eq. (|13|) .as 
well the less stringent condition a > \b 2 \ are violated. 

Fig. 6a shows the difference of x\{t) — x 2 (t — r) Vs t, exhibiting typical feature of on-off 



intermittency 



42 



43| with the off state near the laminar phase and the on state showing 



random burst. In Fig. 6b xi(t) is plotted against x 2 (t — r), where the distribution is scattered 
around the diagonal. To analyze the statistical feature associated with the irregular motion, 



we calculated the distribution of laminar phases A(t) with amplitude less t 



ran a threshold 



42, 



431 ] . where the 



value A = 0.005 as was done in the statistical analysis of intermittency 
power law behavior of mean laminar length is calculated as a function of control parameter. 
A universal asymptotic — | power law distribution is observed in Fig. 6c, which is quite 
typical for on-off intermittency. 

Now, we choose the value of 6 2 = 0.15, below the value of a = 0.16 so that the less strin- 
gent condition a > |6 2 | is satisfied while the general stability criterion Eq. (fT3*|) is violated 
and we carry out the same analysis as above. In Fig. 7a, the difference of x\ (t) — x 2 (t — r) is 
plotted against time t, which is more regular and is much smaller in amplitude but not ex- 
actly zero, thereby implying an approximate anticipatory synchronization Xi(t) m X\(t — r). 
Fig. 7b shows the plot of x\(t) Vs x 2 (t — r), where the distribution is localized entirely 
on the diagonal, but not sharply on it. Earlier we noted that for this case the minimum 
of similarity function S a {r) (Curve 3, inset of Fig. 5) is nearly zero, but not exactly zero. 
The distribution of laminar phase A(t) is plotted in Fig. 7c as for the Fig. 6c. It shows 
a periodic structure in the distribution of laminar phase, where the peaks occur approxi- 
mately at t = nT,n = 1,2, where T is of the order of the period of the lowest periodic 
orbit of the uncoupled system (PJ). It should be remembered that the periodic behavior is 
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associated with the statistical analysis, while the signals remain chaotic. Finally for the 
case 62 = 0.1(< a/1.5), which satisfies the stringent stability criterion ([13)1 . and where the 
similarity function vanishes exactly (Curve 4 in Fig. 5)), exact anticipatory synchronization 
occurs as confirmed in Fig. 4. Thus we find that the transition to exact anticipatory synchro- 
nization precedes a region of approximate anticipatory synchronization from desynchronized 
state as the parameter 62 changes. We have also demonstrated that the emergence of this 
approximate anticipatory synchronization from the desynchronized state is characterized by 
the transition of on-off intermittency to periodic structure in the laminar phase distribution. 

IV. COMPLETE SYNCHRONIZATION FOR r 2 = n 

Complete synchronization follows the anticipatory synchronization as the value of the 
coupling time-delay r 2 equals the feedback time-delay T\, when T2 is increased from a lower 
value. With r 2 = t±, the same stability criterion, Eq. (fT3j) . holds good for this case of 
complete synchronization as well with the same condition 61 = 62 + 03. 

Fig. 8a shows the time trajectory plot of X\(t) and X2(t), exhibiting synchronized evolution 
between them, which is also confirmed by the entirely localized diagonal line of xi(t) Vs 
X2(t) as shown in Fig. 8b. As in the case of anticipatory synchronization, we have found 
that the transition to complete synchronization precedes a region of approximate complete 
synchronization (x\(t) « x^iff) from the desynchronized state as the parameter 62 varies. 
Here also we have identified that the emergence of approximate complete synchronization 
for the case r 2 = T\ is associated with a transition from on-off intermittency to a periodic 
structure in the laminar phase distribution as a function of the parameter 62. In the next 
section we will discuss the existence of lag synchronization for the values of t<i greater than 

n- 

V. LAG SYNCHRONIZATION FOR r 2 > n 

For coupling delay T2 greater than feedback delay Ti, we find that the system (jSJ) exhibits 
exact lag synchronization provided one satisfies the stringent stability criterion (|13|). with 
the lag time equal to the difference between the coupling and feedback delay times. 

Fig. 9a shows the plot of Xi(t) and X2(t) Vs time t, where the response system lags the 
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state of the drive system with constant lag time r = |t2 — t\\. Fig. 9b shows the time-shifted 
plot of Xi(t) and X2 (t + r). However, in the region of less stringent stability condition, 
1.5|6 2 | < a < 1^2 1) approximate lag synchronization occurs as in the cases of anticipatory 
and complete synchronizations. 

We have also calculated the similarity function Si(r) from Eq. (|14j) to characterize the 
lag synchronization. Fig. 10 shows the similarity function Si(t) Vs coupling delay r 2 for four 
different values of &2- Curves 1 and 2 show the similarity function Si(t) for the values of 
62 = 0.18 and 0.16, respectively. The minimum of similarity function Si(t) occurs for values 
of Si(t) > and hence there is a lack of exact lag time between the drive and response 
signals indicating asynchronization. Curve 3 corresponds to the value of 62 — 0.15 (which is 
less than a but greater than a/1.5), where the minimum values of Si(t) is almost zero, but 
not exactly zero (as may be seen in the inset of Fig. 10), so that x\(t) ~ 22 (^ + t). However 
for the value of 63 = 0.1, for which the general condition ()13|) is satisfied, the minimum of 
similarity function becomes exactly zero (Curve 4) indicating that there is an exact time 
shift (Fig. 9) between drive and response signals xi(t) and 22^), respectively, confirming 
the occurrence of lag synchronization. 

We have also confirmed that as in the case of anticipatory synchronization, when the 
parameter 62 varies, the onset of exact lag synchronization is preceded by a region of approx- 
imate lag synchronization, which is characterized by a transition from on-off intermittency 
of the desynchronized state to a periodic structure in the laminar phase distribution. For 
the value of 62 = 0.17 (which violates the stability condition (j!3)l as well as the less stringent 
condition a > I&2I), Fig. 11a shows the difference of x\{t) — X2(t + r) Vs time t, exhibiting a 
typical on-off intermittency. In Fig. lib, Xi(t) is plotted against X2{t + r), where the distri- 
bution is not concentrated along the diagonal. In Fig. 11a, the laminar phase distribution 
A(t) is characterized by an exponential — | power law behavior as shown in the 

Fig. 11c. In order to show that there is a transition from on-off intermittency to periodic 
behavior in the laminar phase distribution corresponding to approximate lag synchroniza- 
tion, we have changed the value of 62 from 0.17 to 0.15, (so that the less stringent condition 
a > 1 6 2 1 is satisfied but not the general condition (jl3|0 . and examined the nature of laminar 
phase distribution A(t). The difference between Xi(t) and X2{t + r) is shown as a function 
of time t in Fig. 12a, where there is only a laminar phase present for the threshold value 
A = 0.002 without any intermittent burst. The corresponding laminar phase distribution 
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A(t) is again characterized by the periodic structure as shown in Fig. 12c. As in the case 
of approximate anticipatory synchronization, here also the peaks occur approximately at 
t = nT, n = 1,2,..., where T is roughly of the order of the period of the lowest periodic orbit 
of the uncoupled system ([TJ. Time-shifted plot x\(t) Vs ^(t + r) is shown in the Fig. 12b, 
where the distribution is concentrated along but not exactly on the diagonal line confirming 
the onset of approximate lag synchronization. As noted previously that for this case the 
minimum of similarity function S\{t) is nearly zero but not exactly zero (Curve 3, inset of 
Fig. 10). Finally for 6 2 — 0.1, which satisfies the general stability criterion (|T3~j) . we have 
exact lag synchronization as demonstrated in Fig. 9 and Fig. 10. Thus we find that as the 
parameter 6 2 varies the transition to exact lag synchronization precedes a region of approxi- 
mate lag synchronization from desynchronized state, where the later is characterized by the 
transition from on-off intermittency to periodic structure in the laminar phase distribution. 



VI. SUMMARY AND CONCLUSION 



In this paper, we have shown the existence of transition from anticipatory synchronization 
to lag synchronization through complete synchronization in a single system of two coupled 
time-delay piecewise linear oscillators with suitable stability condition and with the second 
time-delay r 2 in the coupling as the only control parameter with all the other parameters be- 
ing kept fixed. We have also plotted corresponding similarity functions to characterize both 
the anticipatory and lag synchronization as well as complete synchronization. Further, when 
the parameter 6 2 varies, we find that the transition to exact anticipatory/complete/lag syn- 
chronization precedes a region of approximate anticipatory/complete/lag synchronization 
from desynchronized state, where the region of approximate synchronization is character- 
ized by the transition from on-off intermittency to periodic structure in the laminar phase 
distribution. 

Further, we have observed that in the region where the stringent stability condition (fT3*|) 
is satisfied, the minimum of similarity function S a (r) attains the value zero for all values 
of r 2 < 7~i, indicating that the exact anticipatory synchronization exists for a range of 
coupling delay r 2 below t\. However for approximate anticipatory synchronization (in the 
region 1.5|6 2 | > a > |6 2 |) the minimum of similarity function takes the value S a {r) ~ 0, 
but not exactly zero, for r 2 < T\. Similarly lag synchronization also occurs for a range of 
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delay coupling r 2 above T\. Another interesting aspect is that both the anticipating and 
lag time can be tuned to any desired value by changing the value of coupling delay r 2 . 
Consequently, the kind coupled time-delay systems of the type discussed in this paper have 
considerable physical relevance, particularly for secure communication purposes. We are 
now investigating the existence of similar phenomena in other piecewise linear time-delay 
systems, including the time-delay Chua and Murali-Lakshman-Chua electronic circuits, the 
results of which will be published elsewhere. 
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VII. FIGURE CAPTIONS 

Figl. The schematic form of the piecewise linear function f(x) given by Eq. 
Fig2. The hyperchaotic attractor of the system (0) for the parameter values a = 1.0, 6 = 
1.2 and r = 25.0. 

Fig3. The first ten maximal Lyapunov exponents X max of the scalar time-delay equation 
(JHH) for the parameter values a = 1.0, 6 = 1.2, r G (2, 29). 

Fig4. Exact anticipatory synchronization for the parameter values a = 0.16, b\ = 0.2, b 2 = 
0.1,63 = 0.1, ri = 25.0 and t 2 = 20.0. (a) Time series plot of x\(t) and x 2 (t), (b) Synchro- 
nization manifold between X\{t) and x 2 (t — r), r = r 2 — T\. The response x 2 (t) anticipates 
the drive xi(t) with a time shift of r = 5.0. 

Fig5. Similarity function S a (r) for different values of b 2 , the other system parameters are 
a = 0.16, h = 0.2 and n = 25.0. (Curve I: b 2 = 0.18, 63 = 0.02, Curve 2: b 2 = 0.16, 63 = 
0.04, Curve 3: b 2 = 0.15, 63 = 0.05 and Curve 4: 6 2 = 0.1, 63 = 0.1). 

Fig6. (a) The time series Xi(t) — x 2 (t — r) for b 2 = 0.17 and 63 = 0.03 with all other 
parameters as in Fig. 4 (so that the stability condition is violated for anticipatory synchro- 
nization), (b) Projection of X\(t) Vs x 2 (t — r) and (c) The statistical distribution of laminar 
phase satisfying — | power law scaling. 

Fig7. (a) The time series x\(t) — x 2 (t — r) for b 2 = 0.15 and 63 = 0.05 with all other 
parameters fixed as in Fig. 4 (so that the less stringent condition a > \b 2 \ is satisfied while 
ffT3*|) is violated), (b) Projection of xi(t) Vs x 2 (t — r) and (c) The statistical distribution of 
laminar phase showing a periodic structure. 

Fig8. Exact complete synchronization for the parameter values a = 0.16, 61 = 0.2, b 2 = 
0.1,63 = 0.1, ri = 25.0 and r 2 = 25.0. Here the general stability criterion (|T3|) is satisfied, 
(a) Time series plot of x\(t) and x 2 (t) and (b) Synchronization manifold between Xi(t) and 
x 2 (t). The response x 2 (t) follows identically the drive xi(t) without any time shift. 

Fig9. Exact lag synchronization for the parameter values a = 0.16, 61 = 0.2, b 2 = 0.1, 63 = 
0.1, ri = 25.0 and t 2 = 30.0. Here the general stability criterion p3j1 is satisfied, (a) Time 
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series plot of xi(t) and X2(t), (b) Synchronization manifold between x\(t) and x 2 (t + T). The 
response x 2 (t) lags the drive x\(t) with a time shift of r = 5.0. 

Fig. 10 Similarity function Si(t) for different values of b 2 , the other system parameters 
are a = 0.16, &i = 0.2 and n = 25.0. (Curve 1: b 2 = 0.18, b 3 = 0.02, Curve 2: b 2 = 0.16, & 3 = 
0.04 Curve 3: b 2 = 0.15, b 3 = 0.05 and Curve 4: 6 2 = 0.1,6a = 0.1). 

Figll. (a) The time series x\(t) — x 2 (t + r) for b 2 = 0.17 and 63 = 0.03 with all other 
parameters as in Fig. 9 (so that the stability condition is violated), (b) Projection of xi(t) 
Vs x 2 (t + r) and (c) The statistical distribution of laminar phase satisfying — | power law 
scaling. 

Figl2. (a) The time series Xi(t) — x 2 (t + r) for b 2 = 0.15 and 63 = 0.05 so that the 
less stringent condition a > \b 2 \ is satisfied while f|13j) is violated, (b) Projection of X\{t) Vs 
x 2 (t + r) and (c) The statistical distribution of laminar phase showing periodic structure. 
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